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ABSTRACT 


Many  methods  have  been  proposed  and  tested  in  the 
past  to  obtain  an  optimal  response  for  step  inputs  to 
automatic  feedback  control  systems.  Most  of  these 
methods  have  led  to  the  use  of  sophisticated  control 
devices  ranging  from  small  analog  to  large  digital 
computers.  Here  the  possibility  of  a  simplified 
switching  logic  combined  with  an  open-closed  loop 
servomechanism  is  investigated.  Deadbeat  response  to 
step  inputs  was  the  object  of  this  study  rather  than 
a  time  optimal  response.  Two  types  of  logic  were 
investigated.  A  time  invariant  controller  was  analyzed, 
built  and  tested.  The  system  works  on  the  principle 
of  constant  switching  times  with  the  output  being 
controlled  by  an  open  locp  driving  voltage  which  is 
proportional  to  the  input  step  size.  At  the  completion 
of  the  open  loop  mode  of  operation,  the  system  is 
returned  to  the  normal  closed  loop  mode. 

The  writers  wish  to  express  their  appreciation 
for  the  assistance  and  encouragement  given  them  by 
Dr.  John  Ward,  of  the  U.  S,  Naval  Postgraduate  School, 
in  this  investigation. 
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1 .  Introduction, 

Numerous  schemes  have  been  applied  to  positioning 
feedback  control  mechanisms  to  obtain  optimum  response 
for  step  inputs.  Some  of  the  methods  used  have  been 
(1)  dual  mode  operation  1  using  a  relay  control  computer 
which  approximates  the  optimal  switching  lines,  (2)  use 
of  compensation  networks  to  improve  relay  performance  , 
(3)  application  of  discontinuous  damping  to  a  relay 
servo^,  (4)  conditional  switching  techniques'"1,  and 
others  5,6c  jn  general*  the  basic  idea  behind  these 
various  approaches  is  the  desire  to  obtain  deadbeat 
response  in  minimum  time  for  a  step  input. 

The  time  optimization  requirement  inherently 
requires  a  sophisticated  controller  that  can  compute 
the  proper  switching  times  accurately.  A  question 


*  #  »  #  -X  #  # 

^K,  C.  Matthews,  R.  C.  Boe,  The  application  of 
nonlinear  techniques  to  servomechanisms,  National 
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Nonlinear  Feedback  Control  Systems.  Chapter  7. 
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which  might  he  asked  is:  "If  the  requirements  for 
time  optimization  were  dropped,  retaining  the  provision 
for  deadbeat  response,  would  it  be  possible  to  design  a 
simple  controller  for  a  second  order  system  that  would 
be  of  practical  value?" 

In  this  paper  two  approaches  to  the  development  of 
a  simple  controller  for  a  second  order  system  are 
investigated.  In  both  cases  the  system  operates  in 
two  modes,  linear  and  nonlinear,  which  correspond  to 
closed  and  open  loop  modes  respectively.  In  the  non- 
linear  or  open  loop  mode  a  driving  voltage  is  applied  to 
the  plant  for  a  prescribed  amount  of  time.  The  voltage 
is  then  reversed  until  the  system  reaches  an  output 
which  either  equals  or  approximates  the  input.  At  this 
point  the  controller  returns  the  system  to  the  linear, 
closed  loop,  operation. 

As  previously  stated,  in  both  approaches  investi¬ 
gated,  relay  switching  is  caused  to  occur  on  a  time 
basis.  One  of  these  approaches  sets  the  switching 
times  proportional  to  the  command  signal  while  the  other 
holds  the  switching  times  constant  while  setting  the 
driving  voltage.  The  latter  scheme  is  suggested  by  the 
principle  of  superposition. 

The  system  of  Pig.  1  is  the  basic  block  diagram  for 
both  approaches.  The  step  input  commands  are  represented 
by  R,  the  system  output  position  by  C  and  the  output,  of 
the  open  loop  controller  by  V.  The  plant  is  second  order 

apPao 


and  has  the  transfer  function: 


o) 


K/f 

5(s'Xr) 


The  system  is  to  be  controlled  by  an  open  loop  controller 
which  functions  on  command  signals  only.  Any  load 
perturbations  or  other  similar  disturbances  within  the 
loop  will  not  affect  the  controller. 


Fig,  1  -  Basic  System 

The  assumptions  made  for  the  theoretical  investi¬ 
gations  are  (1)  the  system  is  linear,  (2)  the  system  is 
second  order,  (3)  command  signals  are  to  be  steps  only 
and  (4)  new  commands  are  not  to  be  given  until  the 
previous  command  has  been  completely  executed. 


2.  Time  proportional  to  command  signal. 

With  the  output  voltage  of  an  open  loop  controller 
a  constant  in  magnitude,  there  is  a  specific  relation¬ 
ship  for  the  time  to  reverse  the  polarity  of  the  driving 
voltage  and  the  time  to  return  the  system  to  normal 
linear  operation  in  order  to  force  deadbeat  response. 

By  using  an  approximation  to  this  relationship,  a  theore¬ 
tical  controller  was  derived.  The  system  Incorporating 
this  controller  was  then  investigated  on  the  CDC  digital 
computer. 

For  purposes  of  this  investigation,  the  proposed 
system  is  to  function  as  follows?  At  some  time,  tQ,  a 
step  command  signal,  +  R,  is  applied  to  the  system. 
Application  of  the  command  signal  causes  the  loop  to 
be  opened  and  a  voltage,  +  V,  to  be  applied  to  the  plant 
for  a  period  of  time,  ts.  At  time  tg  the  polarity  of 
the  signal  is  reversed  for  a  period  of  time,  t  .  The 
total  time  of  nonlinear  operation,  tt,  is  defined  by? 

(2)  tt  =  tg  -  tr 

At  time  t^  the  controller  closes  the  loop  returning 
the  plant  to  standard  closed  loop  operation.  If  at 
time  t^,  the  plant  output  has  arrived  at  the  desired 
position  commanded  by  the  input,  then  true  deadbeat 
response  has  been  achieved. 

A  digital  program  was  written  to  compute  the  times 
ts  and  t^  for  various  values  of  inputs.  The  system 
parameters,  K  and  ,  as  well  as  the  applied  voltage,  +V, 
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were  held  constant,.  For  the  system  under  evaluation 
the  following  values  were  arbitrarily  selected? 

K/f  =  1,  l/<f  =  5  and  V  =  100.  Each  combination 
of  tg  and  t^  that  results  in  deadbeat  response  specifies 
a  command  signal  R.  Thus,  from  the  program  data,  a 
plot  of  t0  and  t^  versus  R  was  made  for  05R^TT  radians 
as  shown  in  Fig.  2.  From  these  data,  t  and  t  can  be 
approximated  by  a  straight  line  as  shown  in  the  figure. 
These  straight  line  approximations  of  tg  and  t^  can  be 
evaluated  in  point  slope  form  ass 

(3)  is  =  .075j*\  +  0.05 

(4)  tt  ,092j>?l  +  O.IOS 

Since  tg  and  t^  are  independent  of  the  sign  of  the 
input,  the  magnitude  of  R  is  specified  in  the  above 
equations . 

It  is  immediately  obvious  from  Figc  2  that  the 
straight  line  approximation  is  extremely  poor  for  small 
inputs,  but  this  problem  can  be  eliminated  if  the 
system  is  held  in  closed  loop  for  small  inputs.  Thus, 
in  this  system,  the  minimum  step  size  for  open  loop 
operation  would  be  about  .25  radians. 

Using  the  switching  times  determined  from  equations 
(3)  and  (4),  a  program  was  written  and  appended  to  "Program 
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Fig.  2.  Command  signal  varsue  switching  times 
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Pig.  10.  Transient  response  of 
time  proportional  to  command 
controller.  Input  equals  two 
radians • 


TIME 


7 

Analog"  to  study  the  performance  of  a  system  utilizing 
this  type  controller.  Pour  runs  were  made  with  inputs 
of  1.5»  and  2  radians.  The  resulting  phase 

planes  and  time  response  characteristics  as  plotted  by 
the  computer  are  shown  in  Figures  3  to  10.  By  inspection 
of  Pig.  2  it  is  possible  to  predetermine  whether 
switching  will  occur  early  or  late.  For  examples  if 
the  input  is  ,5  radians 9  Pig,  2  predicts  that  both 
switching  times*  ts  and  tr*  will  be  late  and  the  system 
will  overshoot.  This  is  verified  in  Figures  3  and  4. 

If  the  input  is  1.5  radians*  it  would  be  predicted 
that  switching  would  be  early.  This  is  verified  in 
Figures  7  and  8. 

This  approach  does  offer  some  improvement  in 
time“tO"Steady-state  over  normal  linear  operation  but 
it  is  limited  to  a  comparatively  small  range  of  inputs, 

A  controller  that  could  vary  the  times  of  switching  as 
specified  by  equations  (3)  and  (4)  would  not  be  simple 
to  construct*  which  is  contrary  to  one  of  the  main 
design  objectives  as  set  forth  in  section  1. 

##*##*#* 

^"Program  Analog"  was  written  by  Dr,  J»  R«  Ward 
of  the  U,  S.  Naval  Postgraduate  School,  The  program  as 
used  is  shown  in  Appendix  I. 


ea  |  ^ca 
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Time  invariant  controller 


Assume  that  a  given  linear  second  order  system  is 
controlled  in  such  a  manner  as  to  provide  a  deadbeat 
response  for  a  step  Input.  If  the  output  of  the  linear 
system  is  proportional  to  a  step  input*  then  by  the 

O 

principle  of  superposition  *  the  output  is  doubled 
when  the  input  is  doubled.  Furthermore*  if  the 
acceleration  time*  tc,  *  and  the  deceleration  time  t  , 
are  known  for  a  specified  step  command*  then  it 
will  be  possible  to  obtain  deadbeat  response  for  all 
step  commands  if  the  switching  times  are  held  constant 
while  the  driving  voltage  is  varied  in  proportion  to 
the  step  command. 

To  determine  the  appropriate  switching  times*  it 
is  first  noted  that  if  the  controller  voltage  is  V* 
then  the  system  output  will*  If  the  Initial  conditions 
are  zero,  be» 


’3) 


C(s)  - 


kvA_ 

S2(S  +  /*) 


The  step  voltage  +V  is  to  be  applied  for  a  finite  time, 
tg,  and  the  step  voltage  “V  Is  applied  for  the  finite 


time,  t 


If  at  ts  +  the  output  of  the  system  is 


r*  —  ’'s  t 
identically  equal  to  the  Input*  and  the  output  rate  is 

zero,  then  system  response  is  Indeed  deadbeat.  It  is 


»  #  «  *  #  # 

O 

M.  E.  Van  Valkenburg,  Network  Analysis*  Prentlce~ 
Hall*  Inc.,  pp.  79°80. 
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now  necessary  to  find  the  required  conditions  for  such 
deadbeat  response.. 

The  inverse  transform  of  equation  (5)  is; 

(6)  CM  =  "  0 

To  obtain  the  system  rate  at  time  t  it  is  only  necessary 

to  take  the  derivative  of  equation  (6).  Thus* 

(7)  Cm  s  KV( i  -  £  ^ ) 


The  reversed  voltage  is  now  applied  for  a  period  tr9 

and  considering  the  initial  conditions  that  exist  at 

time  t  9  it  follows  that  the  equations  for  system 
s 

position  and  rate  for  ts  Ik  t  -  t^  are; 

{ 8 )  Ch)  -  VKA  CM  CM  C m/< 

s*Ys  +  y>f)  (s  +  \ /f)  s(s  *  i/f)  S(s  +  y#) 

and; 


(9)  C (*)  s  _  VK/tf  C(tA 

5(5  +  Mr)  ($*Vf) 


or; 


(10) 


O.)  VK 


and; 


(11  ) 


C(h)  ~  VK 


For  deadbeat  response  the  system  rate^  C  » 
must  equal  zero  at  time  t-: ,  Therefore*,  by  setting 
equation  (11)  to  zero*  the  exact  relationship  between  the 


*=  1 7= 


acceleration  time,,  t  „  and  deceleration  time?  t  „  is 


obtained, 

namely? 

(12) 

.-/r  _  £-u« 

£  —  /l-  c 

To  obtain  the  deceleration  time  required  for  any 
acceleration  time,  it  is  only  necessary  to  take  the 
natural  logarithm  of  both  sides  of  equation  (12); 


(12a) 

ir  *  rJ>„(z-e','"f) 

To  obtain  the  system  output  for  deadbeat  response 
in  terms  of  tg  and  tr,  equation  (12)  is  substituted 
into  equation  (10)  which  yields? 

(13)  CM  '  VK(ti-ty) 

Finally,  the  system  output  in  terms  of  the  acceleration 
time  only  can  be  found  by  substituting  equation  (12a) 


Into  (13) 

to  give? 

(14) 

CM  -  VK\t>  -  tfJ.Z-C'  /f\\ 

Thus 

by  using  a  controller  voltage  proportional 

to  R,  and  by  choosing  ts  and  tr  as  calculated  from 
equations  (12a)  and  (14),  deadbeat  response  can  be 
achieved  for  any  step  input  provided  that  the  plant 
Is  adequately  described  by  a  (linear)  second  order 
transfer  function, 

=>  1 8= 


4.  A  design  example  of  the  time  invariant  controller® 

A  logical  sequence  of  manipulations  to  determine 
the  proper  switching  times  for  a  given  linear  second 
order  system  can  he  easily  obtained  from  equations 
(12a)  and  (14).  For  purposes  of  illustration*,  assume 
that  the  system  has  the  following  characteristics? 

K  =  10  Maximum  output  desired  =10 

^=1.0  Saturation,  will  occur  if  KV  is  greater 
than  100® 

Using  equation  (14)  with  the  constants  given  above 
and  selecting  a  maximum  applied  controller  voltage 
of  10  (to  avoid  the  saturation  non-linearity)  a  trans¬ 
cendental  equation  for  the  acceleration  time*,  tg9  is 
obtained? 

(14)  iO  =  tOo\i*  -  &(z -(’**) 1J 

A  value  of  tg  equal  to  ®368  seconds  will  satisfy  the 
equation.  It  is  then  necessary  to  find  the  deceleration 
time  tr*  To  do  this*,  substitu**  time  ts  into  equation  (12a) 

(12a.)  ir  ~  **(  2  ***• 

The  total  time  of  open  loop  operation  is  then  ,6 36  seconds. 

Using  the  above  constants*,  the  system  equations  were 
solved  on  the  digital  computer,  Only  one  input  was 
given  to  the  system  as  other  inputs  would  provide  the 
same  type  of  output  changed  by  a  constant  only.  Fig,  11 
shows  the  computed  phase  plane  response  of  the  system. 
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Test  of  the  time  invariant  controller  applied  to 
a  d-c  motor  plant* 

The  controller  analyzed  in  sections  3  and  4  above, 
is  here  Investigated  experimentally  in  relation  to  a 
DC  motor  plant  which  was  predominantly  second  order 
and  linear  up  to  the  saturation  limit  of  the  driving 
amplifier •  Runs  were  made  with  various  step  inputs  in 
both  the  time” invariant  and  standard  closed  loop  modes* 

The  following  operations  by  the  open  loop  controller 
are  required  to  test  the  theory: 

1)  At  time  ic  a  step  input  is  commanded*  The 
loop  is  opened  and  an  open  loop  driving  voltage  is 
applied  to  the  plant* 

2)  At  time  ts  the  driving  voltage  is  reversed® 

3)  At  time  the  driving  voltage  is  removed 
from  the  plant 9  the  loop  is  closed 9  and  the  system 
returned  to  its  standard  closed  loop  mode  of  operation. 

The  operations  required  were  carried  out  by  a  timing 
device  consisting  of  a  variable  speed  D*  C*  motor 
geared  to  a  shaft  containing  three  “'make  and  break" 
switches.  The  time  of  contact  of  these  switches  rela~ 
tive  to  each  other  was  variable®  that  is®  their  relative 
anglular  positions  on  the  shaft  could  be  varied*  The 
switches  were  connected  to  relay  circuits  which  accomplished 
the  desired  switching*  Thus®  relay  switching  times 

could  be  varied  by  adjusting  the  speed  of  the  D®  C. 
motor  or  by  adjusting  the  relative  angular  position  of 
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the  "make  and  break"  contacts » 

Fig*  12  shows  a  block  diagram  of  the  experimental 
D.  C*  servo  used  for  testing  *  When  a  step  command  is 
received  9  relay  one  switches  to  its  normally  open 
position*  This  is  time  t0»  Relay  two  switches  to  its 
normally  open  position  at  time  tg  to  reverse  the  open 
loop  driving  voltage*  A  third  relay 9  not  shown  in  the 
figure „  is  used  to  return  relays  one  and  two  to  their 
normally  closed  positions  at  time  t+« 


Fig,  12  Co  Servo  with  a  time  invariant  controller 
The  transfer  function  of  the  open  loop  system  was 
evaluated  ass 

07)  iJUJA 

$  (S  +  !/«  9) 

where  the  gain  is  2e94  and  the  system  time  constant  is  1.8 
seconds.  The  system  output  position^,  C9  was  obtained  by 
Integrating  the  output  of  a  tachometer  geared  to  the 
driving  motor.  As  a  result 9  the  system  output  was 
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evaluated  In  volts g  and  consequently 9  the  command  and 
output  signals  could  be  compared  In  volts  rather  than 
radians*  Further  investigation  of  the  open  loop 
system  revealed  that  an  open  loop  command  signal  of 
about  eight  volts  would  cause  saturation  in  the  transistor 
power  amplifier*  A  schematic  of  this  amplifier  is  shown 
in  Fig*  13® 

For  the  purpose  of  steady  state  accuracy 9  of 
primary  importance  in  positioning  sy stems 9  the  closed 
loop  gain  was  set  at  29*4  or  ten  times  greater-  than  the 
open  loop  gain* 

To  simplify  the  construction  of  the  controllers,  it 
was  decided  to  obtain  switching  times  such  that  a  one 
volt  Input*  R*  would  cause  a  one  volt  driving  voltages,  V. 
In  other  words*  the  constant  of  proportionality  derived 
in  equation  (14)  was  made  equal  to  one*  This  means  that 
in  equation  (14)*  V  s  C9  and  tg  can  be  evaluated 
directly  from  that  equation  to  yield  the  desired  value* 


(14) 


-*«/r 


i..  K  e,  -  LBk  ( 2.  -£ 


tt  »  13  sec. 


By  using  the  value  of  tg  equal  to  *98  sec*  In  equation 
(12a) 9  t  is  found  to  be  *54  sec®  The  relay  switching 
times  were  then  set  accordingly* 

An  input  of  three  volts  was  used  to  check  out  the 
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Fig.  13.  Schematic  of  transistor  amplifier. 


switching  time  computed  above*  The  response  was  not 
deadbeat*  The  swl  tailing  times  were  then  adjusted  to 
obtain  deadbeat  response  to  the  three  volt  input*  The 
actual  time  of  voltage  reversal  was  il  0 1 5  sec*  compared 
to  the  computed  ®98  see*  The  total  time  of  operation 
was  1*55  sec*  compared  to  the  computed  time  of  1e52  sec* 
These  differences  were  probably  due  primarily  to  the  fact 
that  the  system  was  only  quasi  linear  and  second.  order 8 
and  the  4#^****  j  Lant  were  determined 

assuming  a  linear  second  order  systems 

A  six  channel  Brush  recorder  was  used  to  record 
the  closed  loop  command  signal 8  the  ©pen  loop  command 
signal^,  the  system  output  and  the  system  rate*  Runs 
were  made  with  the  inputs  varying  from  1*4  volts  to  10 
volts*  A  standard  closed  loop  response  t©  the  same 
step  commands  was  also  made  for  comparison  purposes* 
Figures  14  through  21  show  the  responses  obtained*  Each 
page  contains  a  tlme° invariant  system  response  followed 
immediately  by  a  standard  closed  loop  response  t©  the 
same  signal*  Table  I  is  a  comparison  of  the  different 
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Fig*  l4.  Time  invariant  and  closed  loop  response  of  DC-servo - 

system  with  a  1*94  volt  input*  -  - 
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_ _ _ _ Pig*  15*  Time  invariant  and  closed  loop  response  of  DO  servo 

system  with  a  5*^5  volt  input* 
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'Pig*  16.  Time  invariant  and  closed  lo-’P  response  of  DO  servo 
system  with  a  4.0  volt  input* 
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Fig.  17*  Time  invariant  and  closed  loop  response  of  DCLeervo. 
system  with  a  5* 0  volt  input® 


Fig#  IS#  Time  invariant  and  closed  loop  response  -of  DC  servo  — 
system  with  a  6*p  volt  input* 
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Fig*  .21*  Time  invar iant^and  closed  loop,  response  of  DO  servo., 
system  with  a  1.4  volt  input* 
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The  controller  reversed  the  Input  to  the  motor  at  the 
proper  time  but  the  reversal  Is  not  indicated  on  this 
trace.  It  was  not  possible  to  record  the  Input  to 
the  motor  (amplifier  output)  showing  the  reversal  of 
driving  voltage,,  because  of  grounding  problems,,  A 
ground  at  the  motor  input  would  cause  one  or  both  of 
the  power  transistors  to  burn  out  in  the  power  ampli¬ 
fier.  The  third  trace  is  the  system  output  and  the 
fourth  trace,  the  system  rate.  The  time  ts,  the  time 
at  which  the  polarity  or  sign  of  V  is  reversed^  is 
readily  observed  on  the  last  trace,  system  rate9  as 
there  is  a  definite  change  in  slope  at  that  time® 

Time  on  all  the  traces  reads  from  left  to  right  with  ty 
signified  by  the  commencement  of  the  open  loop  command 
signal  and  tt  by  the  end  of  this  command  and  the  start 
of  the  closed  loop  command  signal. 

It  was  mentioned  earlier  that  an  output  signal  of 
eight  volts  would  saturate  the  power  amplifier  It 
follows,  then,  that  if  ,  the  amplifier  will 

saturate  when  the  driving  voltage  is  reversed.  The 
reason  for  this  is  that  at  the  time  of  voltage  reversal 
a  back  electromotive  force  has  built  up  in  the  motor 
which  approaches  the  driving  voltage.  This  back  emf 
has  the  effect  of  increasing  current  requirement  from 
the  amplifier  when  the  driving  voltage  is  reversed. 
Therefore,  at  the  time  of  switching,  the  system  should 
tend  toward  saturation  if  the  signal  voltage  is  greater 
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than  four  volts*  This  saturation  caused  a  reduction 
in  the  braking  or  deceleration  power  available*  Prom 
this  it  follows  that  the  system  output  would  tend  to 
overshoot.  These  effects  may  be  seen  in  Figures  16, 

17 ,  and  18.  By  analyzing  the  system  rates  at  time  ts 
and  t.j.  for  these  figures  it  can  be  demonstrated  that 
saturation  does  occur  in  the  time  increment  ts  S  t  ‘s  tt« 

At  time  ts,  for  the  three  figures  mentioned  above,  the 
system  rate  increased  proportionately  with  the  command 
signal  as  called  for  by  the  theory.  Yet  at  time 
the  system  rate  was  larger  than  predicted  meaning  that 
the  system  did  not  decelerate  sufficiently  which 
resulted  in  overshoot  at  time  t^.  For  example,  the 
system  rates  at  time  tg  and  for  the  five  volt  input 
(Fig.  17)  are  6.1  and  .95  volts  per  second  respectively. 
For  an  input  of  6.2  volts  (Fig.  18)  one  would  predict 
by  the  principle  of  superposition  that  the  rates  would 
be  increased  by  the  factor  6. 2/5.0  or  would  be  7®5 
and  1.1  volts  per  second.  Actually,  Fig.  18  shows 
the  rates  to  be  7.5  and  1.5  volts  per  second.  The 
rate  at  time  tg  agrees  with  the  theory  but  at  time  t^ 
the  rate  differs  by  .4  volts  per  second.  This  difference 
can  be  accounted  for  by  the  saturation  that  occurs 
when  the  voltage  Is  reversed. 

When  the  system  Input  exceeds  eight  volts,  saturation 
is  predicted  at  t0  and  at  ts  .  One  would  expect,  the re = 
fore,  that  the  system  output  would  never  reach  the 
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desired  point  at  time  t^.  or  in  other  words s  the  system 
will  have  undershoot  at  time  tt»  The  reason  for  this 
phenomenon  is  that  the  driving  voltage  Is  never  large 
enough  to  provide  the  response  called  for  by  the 
theoretical  equations.  The  results  of  this  type  of 
saturation  are  shown  in  Figures  19  and  20.  Fig.  20 
clearly  demonstrates  the  effect  on  system  output. 

Fig.  21  shows  the  response  of  the  system  with  a 
small  input,  1.4  volts.  Here  the  system  output*,  by 
time  tt,  does  not  reach  the  desired  value.  The  primary 
cause  for  this  is  that  the  small  driving  voltage  is  of 
insufficient  magnitude  to  overcome  the  stietion  and 
friction  in  the  system. 

From  these  experimental  tests  it  can  fee  seen  that 
the  time  invariant  switching  scheme  does  provide 
deadbeat  response  over  a  range  of  inputs.  If  deadbeat 
response  is  desired  in  a  positioning  feedback  control 
system,  the  method  of  time  invariant  switching  can 
provide  this  in  a  given  time  increment  over  a  designed 
range  of  inputs. 


6.  Test  of  the  time  invariant  controller  applied  to 

an  amplidyne  driven  motor. 

The  time  invariant  controller  has  been  shown  to 
work  over  a  range  of  Inputs  with  a  quasi  linear , 
second  order  system.  One  might  now  inquire  as  to  its 
usefulness  with  a  system  that  is  quite  nonlinear  in  that 
it  exhibits  hysteresis,  excessive  stietlon,  and  satura- 
tion.  To  investigate  this,  a  system  was  constructed 
that  incorporated  an  amplidyne,  a  i  HP  d“C  shunt  wound 
motor,  and  a  permanent  magnet,  D,  C.  generator.  The 
output  of  the  generator,  or  tachometer,  was  integrated 
to  simulate  the  system  output.  As  in  the  previous 
experiment  the  system  output  and  input  was  in  volts. 

Pig,  22  depicts  the  system  as  set  up  in  the  laboratory 
with  the  relays  performing  the  same  functions  as  out¬ 
lined  in  section  5»  Pig,  23  is  a  picture  of  the  system 
as  it  appeared  in  the  laboratory. 


Pig.  22  Amplidyne  driven  motor  servo  system  with 
time  Invariant  controller. 
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A  plot  of  the  voltage  Input  versus  the  tachometer 
voltage  output  for  the  open  loop  system  appears  on 
Pig.  24.  Notice  that  at  least  two  volts  are  required 
to  start  the  system  in  motion  and  that  at  about  six 
volts  the  system  starts  to  go  into  saturation.  The 
hysteresis  loop  is  rather  broad  which  could  provide 
difficulty  in  system  performance.  Prom  open  and  closed 
loop  tests  of  the  system*  without  the  time  invariant 
switching*  the  open  loop  transfer  function  was  determined 
to  be  approximately? 


(17) 

S(S  t  7) 

Depending  on  the  test  being  conducted*  the  system  gain 
varied  from  100  to  115  and  the  reciprocal  time  constant 
from  6.9  to  7.5  inverse  seconds. 

Using  the  values  of  gain  and  the  time  constant 
from  equation  $7)  and  having  a  unity  relationship  between 
the  open  loop  output  and  open  loop  driving  voltage*  t_ 
was  determined  from  equation  (14)  to  be  about  ,33 
seconds.  Then  t^g  from  equation  (12a)*  was  evaluated 
as  .27  seconds.  The  system  was  rim  with  an  input  step 
of  four  volts.  The  response  was  not  deadbeat  so  the 
times  of  switching  were  adjusted  to  yield  a  near  dead= 
beat  response.  The  timing  circuit  controls  were  not 
fine  enough  to  permit  precise  adjustments  of  the  switching 
times.  The  adjusted  times  were  then  measured  and  found 
to  be  .38  second  for  time  tg  and  .09  second  for  tr. 


“39* 
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Fig.  24.  Hysteresis  loop  of  aaplidyne  driven  servo  systen. 
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The  acceleration  time  was  .05  second  longer  than  the 
computed  value  and  the  deceleration  time  was  .18  second 
shorter  than  the  computed  value.  This  shorter  decelera¬ 
tion  time  is  probably  due  to  the  large  friction  that 
was  present  in  the  motor^ 

A  Brush  recorder  was  used  to  record  the  values  of 
the  closed  loop  command  signal,  the  open  loop  command 
signal,  the  system  output,  the  system  rate,  and  the 
voltage  applied  to  the  motor.  These  recordings  are 
shown  in  Figures  25  through  38*  Following  each  time® 
invariant  run  to  a  given  input,  the  system  was  run 
in  the  normal  closed  loop  mode  with  the  same  quantities 
being  recorded.  A  figure  of  the  closed  loop  recording 
follows  each  time  invariant  figure.  The  maximum  over® 
shoot  and  time~to®steady®state  was  evaluated  for  all 
runs  and  a  comparison  of  the  times  to  steady  state  was 
made.  The  results  of  this  evaluation  are  tabulated  in 
Table  II. 
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2.  Time  to  steady  state  operation  taken  to  be 
when  system  is  within  .1  volts  of  command 
signal. 


3. 


Atgs  is  defined  as  the  difference  between 
the  times  to  steady  state  for  the  open- closed 
loop  and  closed  loop  operation. 


Prom  the  recordings  it  can  be  seen  that  the  range 


from  four  to  five  volts  the  response  was  almost  dead¬ 
beat.  (Figures  25  through  30.)  When  the  input  voltage 


was  in  excess  of  five  volts  the  system  output  started 


to  overshoot  the  desired  value  at  time  t^. 
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Fig.  2 6*  Closed  loop  response  of  amplidyne  driven  motor  servo 
system  with  a  four  volt  input.  . .  . . 
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section  5  with  the  D.  0o  servo 9  this  is  due  to  satura® 
tion  of  the  amplidyne  when  the  driving  voltage  is 
reversed.  This  type  of  response  is  shown  in  Figures  31 
through  34. 

For  step  inputs  less  than  four  volts  system 
performance  deteriorated  rapidly.  This  was  due  almost 
entirely  to  the  excessive  stiction  and  friction  in  the 
system.  It  was  found  that  the  system  did  not  respond  to 
a  command  less  than  two  volts  or  about  30%  of  the 
operating  range.  The  system  response  to  inputs  of  3 
and  3.5  volts  are  shown  in  Figures  35  through  38. 

As  might  be  expected »  the  use  of  an  open^loop 
controller  with  a  far-from= ideal  plant  does  not  lead 
to  acceptable  performance. 
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7.  Conclusions. 

In  the  search  for  a  simple  scheme  to  provide 
deadbeat  response  for  step  inputs  to  a  second  order 
system,  attention  was  concentrated  on  an  open-closed 
loop  controller.  For  the  system  investigated,  switching 
time  in  the  open  loop  mode  were  fixed  and  the  output 
from  the  controller  to  the  plant  was  proportional  to 
the  magnitude  of  the  step  input. 

Using  this  scheme,  any  linear  system  can  be 
controlled  with  deadbeat  response  (Principle  of  super¬ 
position)  .  To  control  the  second  order  plant  considered 
in  this  study  only  one  switching  point  is  required. 

The  experimental  program  was  intended  to  suggest  the 
shortcomings  of  this  approach  when  applied  to  a  relay 
system  with  its  non-linearities  and  higher  order  dynamics. 

The  first  real  system  studied  was  relatively  free 

from  any  non-linearities  and  provided  near  deadbeat 

response  over  approximately  50$  of  the  designed  operating 

range.  A  system  with  rather  pronounced  non-linearities 

was  also  chosen  for  study  and  provided  near  deadbeat 

response  over  less  than  20$  of  the  desired  operating 

range.  It  is  concluded,  then,  that  if  deadbeat  or  near 

deadbeat  response  is  desired,  the  time  invariant  controller 

investigated  in  this  thesis  can  be  used  to  force  this 

response  in  a  linear  system  .  With  a  non-linear  system, 

the  controller  is  of  little  use.  In  any  event,  the  time 

of  response  to  any  size  step  requires  the  same  Interval 

of  time  for  open  loop  operation. 
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APPENDIX  I 


PROGRAM  ANALOG 

A  Control  Data  Corporation  CDC  1604  digital 
computer  is  available  at  the  U*  S.  Naval  Postgraduate 
School.  Dr.  J.  R«  Ward  of  the  Postgraduate  Schools 
Department  of  Electrical  Engineering 9  has  written  a 

O 

FORTRAN  language  program  for  this  computer  which 
simulates  an  analog  computer.  The  program  is  referred 
to  as  Program  Analog*  It  retains  all  the  versatility 
of  the  analog  computer  and  provides  both  a  graphical 
and  numerical  output.  With  this  program  many  of  the 
problems  of  magnitude  and  time  scaling  as  well  as  the 
difficulties  in  simulating  nonlinearities  associated 
with  an  analog  computer  are  avoided.  The  program „  which 
is  reproduced  in  Figures  39  through  41 8  uses  a  modified 
Runge-Kutta  integration  scheme* 

This  program  was  used  for  the  investigations 
carried  out  in  section  2*  Fig.  39  is  the  main  part  of 
the  program.  It  consists  of  the  necessary  Instructions 
for  the  analog  simulation  of  a  given  problem  and  does 
not  require  any  changes  as  the  problem  is  varied.  This 
part  of  the  program  may  be  thought  of  as  the  analog 
computer.  Fig.  40  is  the  FORTRAN  simulation  of  the 

^D.  D,  McCrackens,  A  guide  in  FORTRAN  programmings, 
John  Wiley  &  Sons*,  Inc.8  New  York9  London^  1961. 
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problem  board#  This  part  of  the  program,  must  be 
changed  in  the  same  manner  as  one  would  change  the 
wiring  on  a  problem  board*  Fig*  41 »  the  data  input9 
sets  the  initial  conditions^  potentiometer  settings^ 
problem  time  and  provided  means  by  which  various  output 
data  can  be  obtained*  A  typical  numerical  output 
for  a  Program  Analog  problem  appears  in  Pig*  42* 

Figures  3  through  10  in  the  body  of  the  these  o re 
typical  outputs  for  the  program* 


Fig.  59  •  Program  Analog 


.JOB  MALLEY  (  10  MINUTES  SHOULD  CO) 

PROGRAM  ANALOG 
CALL  EXEC 
ENO 

SUBROUTINE  EXEC 

DIMENSION  X(15C),XCOT(100),C(15C),ITITLE(10),JTITLE(10),KTITLE(40) 

1  ,IT(5),IP(10),IG(10),CT(5),XO(100),PR(10),GR(1C),X1  (900 )  , 

2  Y1  (9  0C )  ,X2 ( 900 ),Y2(900),X3(9C0),Y3(900),X4(9C0),Y4(900), 

3  X5 (90C ) ,Y5 ( 900 ) ,A  !  100  ) 

DO  1  J  = 1 , 1 50 

X  (  J )  =  0. 

1  C  ( J  )  =  0. 

DO  1C01  J=1,10C 

1001  X0( J )  =  0. 

DO  1002  J  =  l,  1C 
I T I TLE ( J )  =  8H 

1002  JTITLE(J)  =  EH 
DO  1003  J=  1 ,40 

1003  KT I TLE ( J )  =  8H 
'  NRC  =  0 

PRINT  200 
200  FORMAT ( 1H1 ) 

2  C 1  FORM ATI / ) 

202  FORMAT!//) 

203  FORMAT! ///) 

READ  102,  (  I  T  ITLE! J  )  ,  J  =  1  ,5) 

102  FORMAT! 10A8) 

READ  101, N, ITEST 
101  FORMAT (  12, A8  ) 

ICHECK  =  8 H  EOUATIC 
IF! ICHECK  -  ITEST)  2,3,2 

2  PRINT  204 

204  FORMAT!  18H  DATA  FORMAT  ERROR) 

STOP  1 

3  READ  101 , NR,  ITEST 
ICHECK  =  8H  RUNS 

IF! ICHECK  -  ITEST )4, 6, 4 

4  ICHECK  =  8H  RUN 

I  Ft  ICHECK  -  ITEST)  7,6,7 

6  I F ( N R  -  9)  3,8,7 

7  PRINT  204 
STOP  2 

8  READ  101 , NT, ITEST 
ICHECK  =  8H  TITLE  C 

IF! ICHECK  -  ITEST)  10,9,10 

9  I F ( NT  -  9)11,11,10 

10  PRINT  204 
STOP  3 

11  DO  12  1=1, NT 

READ  102,  (  JTITLE(J),J=1,  10) 

12  PRINT  205, ( JTITLE! J),J=1 , 1C) 

I F ( N R  -  1)10,13, 14 

13  PRINT  206 
GO  TC  15 

14  PRINT  207, NR 

205  FORMAT ( 5 X ,  1CA3) 

206  FORMAT! 5X.20F1  RUN  IS  CALLED  FOR.) 

207  FORMAT (5X,  II  ,21F  RUNS  ARE  CALLED  FOR.) 

15  PRINT  203 
300  NRC  =  NRC  +  1 

READ  102, ITEST 

I  CHECK  1  =  8HZFRC  CCE 

ICHECK2  =  8HFCLC  CCE 

IF! I  CHECK  1  -  ITEST) 16,  18 , 16 

16  I F (  I CHECK2  -  I  TE ST  )  1  7 , 20 , 1 7 

17  PRINT  204 
STOP  4 

18  DO  19  J=1 , IOC 

19  C(J)  =  0. 

20  RE  AC  10  1  ,NC,  ITEST 


Fig.  59  .  Cont 


21,22,21 


, K=  1  ,  5  ) 


ITEST)  52,32,32 


ICHECK  =  3H  CTEFF. 

I  F (  ICHECK  -  ITEST) 

21  PRIM  2C4 

STOP  5  . 

22  I F (NC) 23, 2d, cZ 

23  DO  24  J=1,MC 
RE  AC  103,  (  IT ( K ) »  CT ( K  ) 

103  FORMAT (5(14, F12. 4)) 

DC  24  Ks 1 , 5 
ITK  =  I T ( K ) 

24  C( ITK)  =  CT ( K ) 

25  PRINT  203, NRC 

208  FORMAT (12H  RLN  NUMBcR 
IRE  AS  FOLLOWS,/) 

K  =  C 

DC  27  J=l,99 
I F (C ( J  )  )26, 27, ct 

26  PRINT  2C9, J  ,  C ( J ) 

209  FORMAT (  1CX2HC (,  12, 4H  )  = 

K  =  K  +  1 

27  CONTINUE  „ 

IF(K)1027,  10^6,1 027 

1026  PRINT  240 

240  FORMAT ( 5X4HNCNE ) 

1027  PRINT  2C2 

READ  102, ITEST  „  r r c 
I  CHECK  1  =  8PZt_RC  ICS 
ICHECK2  =  3HFCLC  ICS 
ICHECK3  =  8HFCLC  FCS 
I F (  I  CHECK  1  -  IT EST ) AC , 28 , 2C 
23  DC  29  J  =  1 , N 

29  XC(J>  =  0 . 

30  F  (  I  CHECK 2  -  ITEST  )3  1^35  ,31  ^ 

31  I  F ( ICHECK3  - 

32  PRINT  204 
STOP  6 

33  DO  3  4  J  =  1 »  N 

34  XC ( J )  =  X ( J ) 

35  READ  101 , N I 
ICHECK  =  3H 
I F (  I  CHECK  - 

'36  PRINT  2C4 
STOP  7 

37  I F  t  N  1)38, 40,33 

38  DO  39  J  =  1  ,NI 
READ  103, ( IT ( K ) 

DC  39  K  = 1 , 5 
ITK  =  I T ( K ) 

39  XC( ITK)  =  CT (K) 

40  PRINT  2  1 C 

210  FORMAT (48H  THE 
K  =  C 

DO  42  J=1,N 
IF ( XC( J ) ) 41  ,  42 , 4  1 

41  PRINT  211 , J,XC( J 

211  FORMAT (  1 0X3HXC ( , I  2 , 4H)  - 

K  =  K  +  1 

42  CCNT INUE 
1  =  ( K  )  1042 , 1  041  ,  1  042 

1041  PRINT  240 
1 C42  PRINT  202 

READ  1C2,iT:_T 
I  CHECK  1  =  BHHCLC  TIv 

ICHECK2  =  8HPEAC  TI-- 
I F ( I  CHECK  1  -  ITcST  43,46 

43  I F I  I CHECK2  -  ITcST)44, 

44  PRINT  204 
STOP  8 

45  READ  104,TC,Tr 


I  1 , / / / , 46H  THE  NON-ZERC  DATA  COEFFICIENTS  A 


E12.5) 


ITEST 
IC  CARC 
ITEST ) 36 


37,  36 
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E12.5) 
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1 04  FORMAT (  10F8. 4 ) 

46  RLAD102 , ITEST 
ICHECK!  =  8HCCMPUTE 
ICHECK2  =  8HF0LC  STE 
ICHECK3  =  8HREAC  STE 

I  F ( I  CHECK  1  -  ITEST  )43, 47 , 48 

47  INDIC1  =  1 

READ  101, KP, ITEST 

I  CHECK  =  8H  IS  CRDE 

I F ( I  CHECK  -  ITEST ) 50 , 59,  50  - 

48  I F (  ICHECK2  -  I T E ST ) 4 9 , 59 , 4 9 

49  I  F (  I  CHECK3  -  I TE ST ) 5 0 , 51  ,  50 

50  PRINT  204 
STOP  9 

51  INDICT  =  0 

READ  104 , (C ! J) , J  =  1 C3 , 109  ) 

I  F  (  C  (  10  3)  -  TO  152,53,52 

52  PRINT  204 
STOP  10 

53  I  F ( C (  104) 154,52,54 

54  NOT  =  1 

I F ( C ( 1 0 5 1  -  TF 1 55 , 59 , 55 

55  I F ( C (106) 156,52,56 

56  NDT  =  2 

I F ( C ( 1 07 1  -TF)57,59,S7 

57  I F ( C ( 1 0 8 1 158,52,58 

58  NDT  =  3 

I  F ( C (  109)  -TF) 52, 59, 52 

59  PRINT  212, TO , TF 

212  FORMAT { 3 1 H  TPE  TIMING  DATA  ARE  AS  FOLLOWS  ,// 

1  5X,  15FINI TIAL  TIME  =  ,E11.5,/ 

2  5X,  151-FINAL  TIME  =  , E 1  1.5) 

IFIINDIC1  -  1)61,60,61 

60  PRINT  213, KP 

213  FORMAT { 5X, 35HTHE  STEP  SIZE  IS  COMPUTED,  BASED  ON  ,/ 

1  5X , 31 FSMALLEST  VARIABLE  CF  CRCER  l.OE  ,12) 

AP  =  10.0**KP* 1 .0E-04/ITF  -  TO 
DO  1060  J=  1 ,  N 
1060  A ( J )  =  AP 

A ( 1 00)  =  TF  -TC 
C { 1 04 )  =  All  CO  1*1. CE-05 
GO  TO  63 

61  NDTT  =  NDT«2 

DO  62  K  =  1 , NDTT, 2 

62  PRINT  214, C( 1C3+K) ,C ( 102+K) ,C( 104+K) 

214  FORMAT! 5X, 15HSTEP  SIZE  =  ,E11.5,10H  FROM  T  =  ,E11.5,3H  TO  T 

1  El  1.5) 

63  PRINT  202 
READ  102, ITEST 
ICHECK  =  8HHCLD  PRI 

I F ( I  CHECK  -  ITEST) 64,74,  64 

64  ICHECK  =  8HREAD  PRI 

IF! ICHECK  -  ITEST)65,1064,65 
1 064  READ  101, NP, ITEST 

ICHECK  =  8H  VARIABL 
IFfICHECK-  ITEST ) 65 ,66, 65 

65  PRINT  204 
STOP  11 

66  IF(NP)67,74, 67 

67  I  F { N  P  -  10)68,68,65 

68  READ  1104  {  I  F  t  K  )  ,K  =  1  ,NP) 

1104  FORMAT (1014) 

READ  101, NT, ITEST 
ICHECK  =  8H  TITLE  C 
IFfICHECK  -  ITEST)65,69, 65 

69  IF(NT)70, 73,70 

70  I F ( N T  -  4)71  ,71 ,65 

71  DO  72  1=1, NT 

K  =  (I  -  1)*10  +1 
KP9  =  K  +  9 
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72  READ  102, (KTITltl J),.J  =  K,KP9) 

73  RE  AC  10  1,  INC  P  R ,  I  TEST 
ICHECK  =  8H  INCREME 
IFIICHECK  -  ITEST >65,74,  65 

74  READ  102, ITEST 
ICHECK  -  8FHCLD  GRA 

I F {  ICHECK  -  ITEST)75,81,  75 

75  ICHECK  =  8HREAD  GRA 

IF (  ICHECK  -  ITEST) 76 ,77,  76 

76  PRINT  204 
STOP  12 

77  READ  101 ,NG, ITEST 
ICHECK  =  8H  GRAPHS 
IFIICHECK  -  ITEST ) 76 ,78, 76 

78  IFING-)  79, 81*79 

79  IFING  -  5)80,80,76 

80  NG2  =  NG*2 

READ  1  104,  {  IGIK  )  ,K=1 ,NG2 ) 

READ  101 ,  INCGR, I  TEST 
ICHECK  =  8H  INCREME 
IFIICHECK  -  ITEST ) 76 ,8 1,  76 

81  PRINT  215 

215  FORMAT! 17H  PRINT  SUMMARY - ,/) 

IF(NP)83,82,S3 

82  PRINT  216 

216  FORMAT ( 5X , 1 1 FNO  PRINTOUT) 

GO  TO  84 

83  PRINT  217, INCPR 

217  FORMAT (5X, I  2 ,29h  INCREMENTS  BETWEEN  PRINTOUTS, 

1  5X ,  25FTHE  VARIABLES  PRINTEC  ARE,/  ) 

PRINT  218,  I  IPI J  )  ,J  =  1  ,NP) 

218  FORMAT! 10X,2FX{ , 13, 1H) ) 

DO  1 C84  J  =  1 , N P 

IF! IPI J) ) 1084, 1 C  85 , 1084 

1084  CONTINUE 
GO  TO  84 

1085  PRINT  1217 

1217  FORMAT  !/,5X,22H*  X ( 0 )  REPRESENTS  TIME  ) 

84  PRINT  202 
PRINT  219 

219  FORMAT! 17H  GRAPH  SUMMARY - ,/) 

IF(NG)86,85, 66 

85  PRINT  220 

220  FORMAT  I 5X,  9FN0  GRAPHS) 

GO  TO  2089 

86  PRINT  221 ,  INCGR 

221  FORMAT! 5X, 12, 26H  INCREMENTS  BETWEEN  POINTS,/  ) 
PRINT  241, I G ! 1 ) , I G ( 2 ) 

241  FORMAT  (  1  OX ,  1 3HGRAPH  A  IS  Xt,I3,8H)  VS.  X ( , I  3 , 
IF! NG  -  2)87,1087, 1087 

1087  PRINT  242,  IG I  3 ) , IG(4) 

242  FORMAT  ( 1  OX ,  1  3HGRAPH  B  IS  Xl,I3,8H)  VS.  X { , I  3 , 
I F ! NG  -  3)87,1088,  1088 

1088  PRINT  243, I G ! 5 ) , IC!6) 

243  FORMAT  { 1  OX,  13HGRAPH  C  IS  X!,I3,8H)  VS.  X  (  ,  I  3 , 
I F ( NG  -  4)87, 1089, 1089 

1089  PRINT  244, I G ( 7 )  ,  I G ( 8 ) 

244  FORMAT  ( 1  OX , 1 3HGRAPH  C  IS  X(,I3,8H)  VS.  X!,I3, 
I F ( N G  -  5)87,1090,1090 

1090  PRINT  245, IGI9) , IG! 10) 

245  FORMAT  ( 1  OX, 1 3HGRAPH  E  IS  X!,I3,8H)  VS.  X  !  ,13, 

87  DO  2087  J=1 ,NG2 

I F ( I G ( J ) ) 2 08 7, 2088, 2087 

2087  CONTINUE 
GO  TO  2089 

2088  PRINT  1217 

2089  PRINT  200 

PRINT  205,  (  ITITLE!  J), J  =  1  ,5) 

PRINT  223, NRC 

223  FORMAT (5X, 10FRUN  NUMBER,  I  2 , / // ) 


1H)  ) 
1  H  )  ) 
1H)  ) 
1H  )  ) 
1H)  ) 
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IF(NP)88,90, 88 

88  I  F (NT)89,90, 89 

89  PRINT  224, (KTITLEf J) ,J=1 ,KP9 ) 

224  FORMAT ( 2  X , 9 ( A  8 , 4X)  ,A3) 

PRINT  201 

90  T  =  TO 

DT  =  C(  104  ) 

DO  91  J=1 ,N 

91  X ( J )  *  XO ( J ) 

LINES  =  0 
NCPTS  =  0 
NUMPTS  =  0 

301  IF(NP)302, 312,302 

302  IF(N0PTS)303,308,303 

303  I  F ( XMCDF(N0PTS»50*INCPR)  1304,306,304 

304  I F ( XMCDF ( NOP  TS ,  10* INC PR)  130  5,307,30  5 

305  I F ( XMODF ( NOPTS  t  INCPR)  1312,308,312 

306  PRINT  200 

307  PRINT  201 

308  CALL  DERIVI T,X, XDOT,C) 

LINES  =  LINES  +  1 

00  311  J  =  1 ,NP 
I  F (  I  P ( J 1 1310  ,309,310 

309  P R ( J  1  =  T 
GO  TO  311 

310  IPJ  =  I P ( J  1 

P  R  (  J 1  =  XI IPJ) 

311  CONTINUE 

PRINT  225, ( PR ( J ) , J= 1 ,NP) 

225  FORMAT (  1  0 (  IX, El  1.5)1 

312  IFING1313, 31 8,313 

313  I F ( XMCDF (NOP TS,  I  NC  GR  )  )  3 1  8 , 3  1  4 . 3  1  8 

314  I F( XMODF (NOPTS, IMCPR) 1 1315, 1314, 1315 
1315  CALL  DERIV(T,X,XDOT,C) 

1314  DO  317  J=1 , NG2 

I F ( I G ( J 1 1316,315,316 

315  GR ( J 1  =  T 
GO  TO  317 

316  IGJ  =  I G ( J 1 

GR ( J  1  =  X(  IGJ) 

317  CONTINUE 

NUMPTS  =  NUMPTS  +  1 
Y 1 (NUMPTS)  =  GR ( 1 1 
XI  (NUMPTS)  =  GR( 2) 

Y2( NUMPTS)  =  GR ( 3 ) 

X2 ( N UMP  TS ) =  GR ( 4 ) 

Y  3 ( NUMPTS)  =  GR ( 5 ) 

X  3 ( NUMPTS)  =  GR ( 6 ) 

Y4 (NUMPTS)  =  G  R ( 7 ) 

X4 ( NUMP  T  S )  =  GR ( 8 ) 

Y5( NUMPTS)  =  GR ( 9 ) 

X5 ( NUMP  TS )  =  GR(  10) 

318  NOPTS  =  NOPTS  +  1 

I  F ( L  I N E S  -  250)1319,  1318,131  8 

1318  PRINT  1216 

1216  FORMAT  (//24F  STOP  AT  250  PRINT  LINES) 
GO  TO  341 

1319  I F ( T  -  1. £  +  04)320, 319, 31  9 

319  PRINT  226 

226  FORMAT ( / / 1 9H  STCP  AT  T  =  10,000) 

GO  TO  341 

320  IFtNOPTS  -  1C000  )322,321,321 

321  PRINT  227 

227  FORMAT ( //26H  STCP  AT  10,000  INCREMENTS) 
GO  TO  341 

322  I F ( NUMPTS  -  900)324,323,223 

323  PRINT  228 

228  FORMAT ( //25H  STCP  AT  900  GRAPH  POINTS) 
GO  TO  341 

324  I F ( T  -  TF) 326,325, 325 
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325  PRINT  229 

229  FORMAT ( / /26H  NORMAL  STOP  AT  FINAL  TIME) 

GO  TO  341 

326  DO  328  J=1 , 1 50 

I F ( A8SF (X ( J)  )  -  1. £  +  04)328,327, 327 

327  PRINT  230, J 

230  F0RMAT(//1 lH  STOP  AT  X(,I3,1CH)  =  1C, COO) 

GO  TO  341 

328  CONTINUE 

330  I F ( I  NO  I C 1  >331. 332*  1 331 

331  STOP  13 

332  CALL  RKUTTA  (N,T,DT,X,C) 

I  F ( NOT  -  1  )334, 333,334 

333  DT  =  C ( 1 04 ) 

GO  TO  340 

334  I F ( NOT  -  2)  337,335,337 

335  I  F ( T  -  C (  105)  )333, 336,336 

336  DT  =  C ( 106 ) 

GO  TO  340 

337  I  F ( T  -  Cl  105) 1333,338,338 

338  I F ( T  -  C( 107)  )35o,339,339 

339  DT  =  C ( 108) 

340  T  =  T+DT 
GO  TC  301 

1331  CALL  RKUTTA2  ( N  ,  T , CT , X , A , C ) 

GO  TO  301 

341  I F(NG) 1341,356,1341 

1341  GO  TO  ( 342,343,344,345,346,347,348,349,350) ,NRC 

342  ITITLE16)  =  £H  RUN  1 
GO  TO  351 

343  ITITLE16)  =  8H  RUN  2 
GO  TO  351 

344  ITITLE16)  =  8H  RUN  3 
GO  TC  351 

345  ITITLE16)  =  8H  RUN  4 
GO  TC  351 

346  ITITLEl  6)  =  £H  RUN  5 
GO  TO  351 

347  ITITLE(o)  =  8H  RUN  6 
GO  TO  351 

348  ITITLE16)  =  £H  RUN  7 
GO  TO  351 

349  ITITLE16)  =  8H  RUN  8 
GO  TO  351 

350  ITITLEl 6 )  =  8H  RUN  9 

-  35 1  LABEL  =  4H 

MCDCLRV  =  0 
SFX  =  0. 

SFY  =  0. 

MINCFFX  =  0 
MINOFFY  =  0 
LABELNO  =  11 
MODE  =  0 

ITITLEl 7 )  =  8H  GRAPH  A 

CALL  GRAPH2  (  NUN PT S , X 1 , Y 1 , 8 , MCDCURV, L ABEL , I T I T L E ,  SFX ,  SF Y, 
1  MINCFFX, MIN0FFY,LA8ELNC»MCDE) 

I F ( N G  -  2)356,352,352 

352  ITITLE17)  =  EH  GRAPH  B 

CALL  GRAPH2  { NUM P T S , X2 , Y 2 , 8 , MODCURV, L ABEL , I T I TLE , SFX ,  SF Y, 
1  MINOFFX, MINOFFY, LABEL NC, MODE) 

I F ( NG  -  3)356,353,355 

353  ITITLEl 7)  =  EH  GRAPH  G 

CALL  GRAPH2  ( NUM PT S , X 3 , Y 3 , 8 , MCDCURV , L ABEL , I T IT Lc , SF X , S F Y, 
1  MINOFFX, MINCFFY,LABELNC,MCDE) 

IF(NG  -  4) 356,354,354 

354  ITITLE17)  =  cH  GRAPH  D 

CALL  GRAPH2  ( NUMPT S, X4 , Y4 , 8 , MCDCURV, L ABEL , I T I TLE  ,  SFX,  SF Y, 
1  MINOFFX, MIN CFFY,LABELNC, MODE) 

I  F ( NG  -  5)356,355,355 

355  I T I T L E ( 7 )  =  EH  GRAPH  E 
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CALL  GRAPH2  (NUMPTS, X5 , Y 5 , 8 , MODC URV , LABEL,  [TITLE,  SFX,  SPY, 

1  MINQFFX,MINCFFY ,  LABELNC , MODE ) 

356  PRINT  200 
I F ( NRC  -  NR) 300,357, 357 

357  RETURN 
END 

,  SUBROUTINE  RKUTTA  (N,T,DT,X,C) 

DIMENSION  X(  150) , AK{4, IOC) ,X COT (100) , XC(  100)  ,CT( 4 ) , C (150 
CT ( 1  )  =  0.0 

C T ( 2  )  =  0.5 
CT ( 3  )  =  0.5 
C  T { 4 )  =  1.0 
DO  4  1=1,4 
TC  =  T  +  CT ( I ) *0T 
DC  2  J= 1 , N 

2  XC(J)  =  X ( J )  +  CT(  I)*AK(  1-1, J) 

CALL  DERIV  (TC,XC,XDOT,C) 

DO  4  J= 1 , N 

4  AK (  I  , J )  =  DT  «XDCT ( J ) 

DO  3  J  =  1 , N 

3  X(J)  =  X(J)  +  (AK( 1,J)+2.*AK(2, J)+2.*AK(3,J) +AK(4,J )  )*0.1  6666  66667 
RETURN 
END 

SUBROUTINE  RKUTTA2  ( N , T, CT , X , A , C ) 

THIS  SUBROUTINE  IS  BASEC  ON  THE  METHOD  DESCRIBED  IN  COMM.  OF 
THE  A.C.M.,  JUNE  I960  (PP.  355  -360).  IT  CHECKS  TFE 
TRUNCATION  ERROR  AT  EACH  STEP  OF  THE  INTEGRATION  AND  ACJUSTS 
THE  STEP  SIZE  ACCORDINGLY.  THE  ACTUAL  RUNGE-KUTTA  INTEGRATION 
IS  PERFORMED  EY  SUBROUTINE  RKUTTA,  WHICH  IS  AVAILABLE  IN  THE 
USNPGS  COMPUTING  CENTER.  THE  ARGUEMENTS  ARE, 

N  =  NUMBER  OF  (FIRST  ORCER)  EQUATIONS.  MAXIMUM  N  =99. 

T  =  TIME  AT  START  OF  INTEGRATION  STEP  (UPDATED  BY  THIS 
SUBROUTINE  AFTER  THE  COMPLETION  OF  EACH  STEP). 

DT  =  STEP  SIZE  AT  START  OF  EACH  STEP  (ALSO  UPDATED). 

X ( I )  =  THE  N  DEPENDENT  VARIABLES  (ALSO  UPDATED). 

A(I)  =  THE  SPECIFIED  ALLOWABLE  ERROR  PER  UNIT  TIME  FOR  EACH 
OF  THE  DEPENDENT  VARIABLES. 

A (  100)  IS  USED  TC  ENTER  THE  TOTAL  TIME. 

NOTE  THAT  IF  NECESSARY  DT  IS  REDUCED  FROM  THE  VALUE  STATED 
IN  THE  ARGLEMENT  UNTIL  THE  SPECIFIED  ACCURACY  HAS  BEEN 
ACHI EVEC. 

DIMENSION  X (  150  )  ,XS(  ICO)  , X22 ( 1  CO  )  , X2 (  1  00 ) , A (  1 C 0  )  , C (  1  50 ) 

TS  =  T 
H  =  2.0*DT 
DO  1  I  =  1,N 
X2( I  )  =  X  (  I  ) 

1  X22  (  I  )  =  X  (  I  ) 

CALL  RKUTTA  ( N , T S , H , X2 2 , C ) 

H  =  DT 

2  CALL  RKUTTA  (N,TS,H,X2»C) 

DO  3  1=1 , N 

3  XS (  I  )  =  X2 ( I  ) 

TS  =  TS  +  H 

CALL  RKUTTA  ( N , T S , H , X S ,C ) 

U2  =  0.03 
DO  6  1=1, N 

E  2 1  =  { X22 ( I )  — X  S ( I ) )*0.  06666666667 

E  2 1 R  =  DIMF ( ABSF ( E21 ) , AB SF ( ( ABSF ( XS ( I ) ) +ABSF ( X22 ( I ) ) 

1  -1 .99*ABSF (X ( I ) ) ) *1 .0E-C9 ) ) 

C  THIS  CONDITION  PREVENTS  RCUNC-OFF  ERROR  FROM  TAKING  CONTROL. 

U  =  E21R  /  ( A ( I)*6.0*H) 

IF  (U  -  U2 ) 6 , 6 , 5 

5  U2  =  U 

6  X S (  I  )  =  XS ( I  )  -  E21 
IF  ( U2  -  1.0)11 ,7,7 

7  IF  (H  -  A ( IOC )* 1 .0E-09 )8 ,8 ,9 

8  DT  =  A( 100)*1.0E-09 

C  THIS  SETS  THE  MINIMUM  STEP  SIZE  TO  1.0E-09  TIMES  THE  TCTAL  TIME. 

GO  TC  16 

9  DO  10  I =  1 ,N 
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X22 (  I)  =  X 2 {  I  ) 

10  X  2  t  I  )  =  X  (  I  ) 

TS  =  T 

H  =  C.5*H 
GO  TO  2 

THIS  RECYCLES  THE  INTEGRATION  IF  THE  TRUNCATION  ERROR  IS  EXCESSIVE. 

11  IF  ( U2  -  0.031 ) 12, 13, 13 

12  DT  =  2.0*H 
GO  TO  14 

13  DT  =  SQRTFt SCRTF (0.5/U2)  )*b 

14  IF  (DT  -  A( ICO) <0.001 ) 16 ,16,15 

15  DT  =  A (  100) *C.CC  1 

THIS  SETS  THE  MAXIMUM  STEP  SIZE  TO  1.0E-03  TIMES  THE  TOTAL  TIME. 

16  T  =  T  +  2 . 0*  H 
C  (  1C4 )  =  DT 
DO  17  I  = 1  ,  N 

17  X ( I )  =  XS  (  I  ) 

THIS  UPDATES  T  AND  X ( I ) .  DT  IS  UPDATED  BY  STATEMENT  8,12,13  OR  15. 

RETURN 

END 

FUNCTION  RELAY  (R,DZONE,V) 

I F ( A BSF ( R  )  -  DZCNE)3,3,4 

3  RELAY  =  0.0 
RETURN 

4  RELAY  =  S I GNF ( V , R ) 

RETURN 

END 

SUBROUTINE  DERIV  (T,X,XDCT,C) 

DIMENSION  X (  150)  ,XDOT( 10  0),C( 150) 


noon  ooo 


Fig.  40.  Subroutine  Deriv  for  Program  Analog 


COMMENTS 

SUBROUTINE  FCR  CPEN  -  CLOSED  CYCLE  SYSTEM  USING  A  TIMING  C  K  T  WHERE 
TIME  TO  SWITCH  AND  TIME  TC  TURN  OFF  IS  PROP  TC  THE  COMMAND  SIGNAL 
AND  THE  RELAY  OUTPUT  IS  A  CONSTANT  IN  THE  NONLINEAR  REGION 
R  =  C  (  1  ) 

ERROR  =  R  -  X ( 1 ) 

IF(SWl)  599,  599,  600 


599 

IF  { ABS F ( ERROR  -  C(2)  )  ) 

603, 

600  , 

600 

600 

SIGN  = 

ABSF ( ERROR ) /ERROR 

SW1  = 

1.0 

SW2  = 

O.C 

IF  (  T 

-  TS  -  TIME)  601 , 

602, 

602 

601 

DRIVE 

=  C ( 6 )  *  SW1 

GO  TO 

604 

602 

DRIVE 

=  -C { 6 )*  SW1 

I F ( ABSF ( XDOT ( 1 ) )  -  .05) 

603, 

603, 

604 

603 

SW1  = 

0.0 

S  W2  = 

1.0 

TIME  = 

X(0) 

604 

POWER 

=  C ( 3 )  *  ERROR  *  SW2  + 

DRIVE 

*  SWl 

XDOT (1 )  =  X ( 2 ) 

XDOT ( 2 )  =  POWER  *  C(4)  -  X { 2 )  *  C ( 5 ) 

TS  =  0. 075  *  ABS  F ( R )  +  0 .05 
TO  =  0.092  *  ABSF(R)  +  0.105 
X ( 3 )  =  SW1 

X  ( 4 )  =  SW2 
X  (  5 )  =  DRIVE 
X ( 6 )  =  POWER 
X(lOl)  =  ERROR 
X( 102)  =  —XDOT (  1  ) 

RETURN 


0 MM ENTS  THIS  IS  BLUE  DECK  ONE 

'  X  (  1  ) =SYSTEM  OUTPUT ,  X(2)=SYSTEM  RATE,  C(l)=  INPUT, 
-  OF  LINEAR  ZONE,  C(3)=SYSTEM  GAIN,  C(4)=PLANT  GAIN, 
C ( 6 ) =MA GN I TUCE  OF  DRIVING  VOLTAGE  IN  NONLIN  REGION 
END 


C ( 2  )  =  1/2WIDTH 
C { 5  )  =  1 /TIME  C 


Fig*  4l-  Data  input  for  Program  Analog 


Kill tKMAN  NALLE Y 
02  EQUATIONS 
04  RUNS 

01  TITLE  CARDS 
OPEN-CLOSED  LOOP 
ZERO  COEFFS 
06  COEFF.  CARCS 


IHES  I  S 


SERVC  SYSTEM  WITH  SWITCHING  =  F ( T I  ME , 


1 

2 

3 

.  4 

5 

6 

ZERO 
00  IC 
READ 
0.0 


0.5 
0.25 
40.00 
1  .0 
5.0 
100.0 
ICS 
CARDS 
TIME  DATA 
1 . 5 


COMPUTE  STEP  SIZE 

-1  IS  ORDER  OF  SMALLEST  VARIABLE 
READ  PRINT  CATA 
07  VARIABLES  PRINTED 

0  1  2  3  4  5  6  7 

02  TITLE  CARDS 

TIME  SYSTEM  SYSTEM  SWITCH  SWITCH 
OUTPUT  RATE  ONE  TWO 

20  INCREMENTS  BETWEEN  PRINTS 
READ  GRAPH  DATA 
02  GRAPHS 

1  0  102  101 

01  INCREMENTS  BETWEEN  POINTS 
HOLD  COEFFS 
01  COEFF.  CARDS 
1  1.0 
ZERO  ICS 
00  IC  CARDS 
HOLD  TIME  DATA 
HOLD  STEP  SIZE 
HOLD  PRINT  DATA 
HOLD  GRAPH  DATA 
HOLD  COEFFS 
01  COEFF.  CARDS 
1  1.5 

ZERO  ICS 
00  IC  CARCS 
HOLD  TIME  DATA 
HOLD  STEP  SIZE 
HOLD  PRINT  CATA 
HOLD  GRAPH  DATA 
HOLD  COEFFS 
01  COEFF.  CARCS 
1  2  .0 
ZERO  ICS 
00  IC  CARDS 
HOLD  TIME  CATA 
HOLD  STEP  SIZE 
HOLD  PRINT  DATA 
HOLD  GRAPH  DATA 


CRI  VE 


PCWER 
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INPUT ) 


Fig.  42.  Typical  output  for  Program  Analog 


OPEN-CLOSED  LOOP  SERVO  SYSTEM  WITH  SWITCHING  =  F ( T I  ME »  INPUT) 
4  RUNS  ARE  CALLED  FOR. 


RUN  NUMBER  1 

THE  NON-ZERO  DATA  COEFFICIENTS  ARE  AS  FOLLOWS 


C( 

1 ) 

= 

. 50000  E+OD 

C( 

2) 

= 

.25000E+00 

C( 

3) 

S 

.4  0000E  +  02 

c< 

4) 

— 

. 1 0000E+01 

C( 

5) 

= 

. 50000  E  +  01 

C( 

6) 

= 

. I0000E+03 

THE  NON-ZERO  INITITAL  CONDITIONS  ARE  AS  FOLLOWS 
NONE 


THE  TIMING  DATA  ARE  AS  FOLLOWS 

INITIAL  TIME  =  .OOOOOE+DO 
FINAL  TIME  =  .15000E+D1 

THE  STEP  SIZE  IS  COMPUTED,  BASED  ON 
SMALLEST  VARIABLE  OF  ORDER  l.OE-1 


PRINT  SUMMARY - 

5  INCREMENTS  BETWEEN  PRINTOUTS 
THE  VARIABLES  PRINTED  ARE 


X( 

0) 

X( 

1 ) 

X( 

2) 

X( 

3) 

X( 

4  ) 

X( 

5) 

X( 

6) 

*  X(O)  REPRESENTS  TIME 

GRAPH  SUMMARY - 

NO  GRAPHS 
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KITTERMAN  MALLEY  THESIS 
RUN  NUMBER  1 


TIME 


.OOOOOE+OO 
. 93000E- 03 
.128108-01 
.278 1 0E-01 
. 42  8 1 OE-0 1 
.578 1 OE-O 1 
.7281 OE-O I 
.8781 0E-01 
.  10281E  +  00 
.  1  1781E+00 


.  1 3281 E+00 
.  1  4500E  +  00 
. 14500E+00 
. 1 4500E+00 
.  14500E+00 
.  14502E  +  00 
. 14579E+00 
.  1 57 1 5E  +  00 
. 1 72 1 5E+00 
.■18715E  +  00 


.202 1 5E  +  00 
.2171 5E+00 
. 232 1 5E+00 
.2471 5E+00 
.262158+00 
.2771 5E  +  00 
. 292 1 5E+00 
. 307 1 5E  +  00 
.3221 5E+00 
.3371 5E  +  0C 


.35215E-S-00 
.3671 5E+00 
.  382 1 5E  +  00 
. 397 1 5E+  00 
.4121 5E  +  00 
.427 1 5E+00 
.4421 5E+00 
. 4571 5E+00 
.4721 5E+00 
.4871 5E+00 


.  502 1 5E  +  00 
.  5  1 7 1 5E  +  00 
.53215E+00 
.54715E+00 
.  562 1 5E  +  00 
.5771 5E+00 
.  5921 5E  +  00 
.60715E+00 
.6221 5E  +  00 
.6371 5E+00 


SYSTEM 

OUTPUT 


.OOOOOE+OO 
-.43173E-04 
-. 80324E-02 
-. 36938E-01 
-. 85432E-01 
-. 1521 OE+OO 
-. 23563E+00 
-. 33479E+00 
-. 44847E+00 
-. 57562E+00 


-. 71 525E+00 
-. 83727E+00 
- . 8  3730E  +00 
83730E+00 
- . 8373 1 E+00 
-.83755E+00 
-. 84538E+00 
- . 95485E  +  00 
-.10790E+01 
-. 1 1805E+01 


12601E+01 
-. 13187E+01 
-. 13574E+01 
-.  1  3771 E  +  0 1 
-. 1 3791 E+01 
-. 13648E+01 
-. 1 3353E+0 1 
1 2920E+01 
-. 12363E+01 
-.  1  1696E+01 


-. 10933E+01 
-. 10087E+01 
-.91707E+00 
-. 81 983E+00 
-.71817E+00 
-.61331 E+00 
-.50638E+00 
-. 39845E+00 
- . 2  9054  E+  00 
-. 18358E+00 


-.784228-01 
.241498-01 
. 1 2343E+00 
. 2 1 880E+00 
. 30972E+00 
. 3957 1 E+00 
. 476  40  E  +  00 
. 551 46E+00 
. 62065E+00 
.683808+00 


SYST  EM 
RATE 


.OOOOOE+OO 
-.927  84E-0 1 
-. 12408E+01 
- . 259  63E  +0 1 
-.38538E+01 
- . 502  05  E  +  0 1 
-.6 1029E+01 
-. 7 1 0  70E  +  01 
-.803  86E  +  01 
-. 89029E+0 1 


- . 970  48E  +  0 1 
-.  10313E+02 
-.103 1 4E+02 
-.1031 4E+02 
-.103 1 3E+02 
-.10311 E+02 
-.  10231  E+C2 
- . 904  69E  +0 1 
-.751  48E+01 
-. 602  82  E+0 1 


-. 45969E+01 
-.32291 E+01 
-. 19322E+01 
-.712  32E  +  00 
.425  63E+00 
.14777E+01 
.244098+01 
.331  34E  +  01 
.  409  4  1  E  +  0 1 
. 473  29E  +  0 1 


. 533  05E  +  01 
.583  84E  +  C 1 
.  6  30  85E  +  0 1 
. 664  36E  +  01 
.689  69E+0 1 
.707  22  E+0 1 
.717  35E+0 1 
.720  53E  +  01 
.71725E+01 
.707  98E  +  01 


. 693  25E  +  01 
. 673  58E  +  0 1 
.649  49 E  +  0 1 
.621 52E+01 
.590 18E+01 
. 553  99E  +  0 1 
.519 46 E+01 
.481 07E+01 
.441  28E  +  01 
.400  54E  +  0 1 


SWITCH 

ONE 


. OOOOOE+OO 
. 1 OOOOE+O 1 
. 1 OOOOE+O 1 
. 10000E+01 
. 10000E+01 
. 1 OOOOE+O 1 
. 10000E+01 
.  1  OOOOE  +  O 1 
. 10000E+01 
. 10000E+01 


. 1 OOOOE+O 1 
. 1 0000E+01 
. 1 OOCO  E+0 1 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 


.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
. OOOCOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 


.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
. OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 


.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
•.OOOOOE  +  OO 
.OOOOOE+OO 
.OOOOOE+OO 
. OOOOOE+OO 


SWITCH 

TWO 


.  10000E  +  01 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
. 00000c  -  90 
. OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 


.OOOOOE+OO 
.OOOOOE+OO 
.OOOOOE+OO 
.  1 0000  E  +  0 1 
. 10000E+01 
. 1 OOCOE+O 1 
. 1 OOOOE+O 1 
. 1 OOOOE+O 1 
. 10000E+01 
. 1 OOOOE+O 1 


. 10000E+01 
. 10000E+01 
. 1 OOOOE+O 1 
. 10000E+01 
.  10000E  +  01 
. 1 OOOOE+O 1 
.  1  0  000  E+  0 1 
. 1 0000  E+0 1 
. 10000E+01 
. 1 OOOOE+O 1 


. 10000E+01 
. 10000E+01 
. 1 OOOOE+O 1 
. 10000E+01 
.  1 0000  E  +  0 1 
. 10000E+01 
. 1 0000E+01 
. 1 OOOOE+O 1 
. 10000E+01 
. 1 OOOOE+O 1 


. 10000E+01 
. 10000E+01 
. 1 OOOOE+Ol 
.  1  OOOOE  +  Ol 
. 1 OOOOE+O 1 
. 10000E+01 
. 10000E+01 
. 1 OOOOE+O 1 
. 1 OOOOE+Ol 
. 10000E+01 


DR  I  VE 


.  1  OOOOE  +  O  3 
. 10000E+03 
.  1 OOOOE  +  O  3 
. 1 0000E+03 
. 10000E+03 
.  10000E  +  03 
. 1 00005+03 
. 1 OOOOE  +  O  3 
, 1 0000E+03 
.  1  0000E  +  03 


.  1 0000  E  +  0  3 
.  1  0000E  +  03 
.  10000E  +  03 
.  1 0000  E+0  3 
. 1 OOOOE+03 
. 1 0000E+03 
. 1 OOOOE+03 
. 1 0000E+03 
. 1 0000E+03 
. 1 0000E+03 


. 1 0000E+03 
. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE  +  O  3 
.  1 OOOOE  +  O  3 
. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE  +  O  3 
. 1 OOOOE  +  O  3 
. 1 0000E+03 


. 1 0000E+03 
. 1 0000E+03 
.  1  OOOOE  +  Cr 
. 1 0000E+03 
.  1  OOOOE+O  3 
. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE+O  3 
.  1  0000E+03 
.  1  0000E+03 


. 1 0000E+03 
. 1 0000E+03 
. 1 0000E+03 
. 1 OOOOE  +  O  3 
. 1 0000E+03 
. 1 0000E+03 
.  1 OOOOE  +  O  3 
. 1 OOOOE+O  3 
. 1 0000E+03 
. 1 OOOOE+03 


i. 
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